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ABSTRACT 

We have used the Hubble Space Telescope's Advanced Camera for Surveys dFord et al . 2003) to mea- 



sure the cumulative mass density in morphologically-selected early-type galaxies over the redshift range 
0.8 < z < 1.7. Our imaging data set covers four well-separated sight-lines, and is roughly intermediate 
(in terms of both depth and area) between the GOODS/GEMS imaging data, and the images obtained in 
the Hubble Deep Field campaigns. Our images contain 144 galaxies with ultra-deep spectroscopy ob- 
tained as part of the Gemini Deep Deep Survey. These images have been analyzed using a new purpose- 
written morphological analysis code which improves the reliability of morphological classifications by 
adopting a 'quasi-Petrosian' image thresholding technique. We find that at z ~ 1 about 80% of the stars 
living in the most massive galaxies reside in early-type systems. This fraction is similar to that seen in the 
local Universe. However, we detect very rapid evolution in this fraction over the range 0.8 < z < 1.7, 
suggesting that over this redshift range the strong morphology-mass relationship seen in the nearby Uni- 
verse is beginning to fall into place. By comparing our images to published spectroscopic classifications, 
we show that little ambiguity exists in connecting spectral classes to morphological classes for spectro- 
scopically quiescent systems. However, the mass density function of early-type galaxies is evolving more 
rapidly than that of spectroscopically quiescent systems, which we take as further evidence that we are 
witnessing the formation of massive early-type galaxies over the 0.8 < z < 1.7 redshift range. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 

The study of galaxy formation and evolution is one 
of the most active interfaces between observation and 
theory in contemporary astrophysics. Sophisticated 
numerical simulations have elucidated the key role 
of d ark matter in driving the formation of structure 
(e.g. Springel et al. I2005I) . while the most recent gen- 
eration of simulations incorporate gas dynamics and 
highlight the role of fe edback (e.g. Croton et al.| [2006; 
Governato et al l2006h . A number of recent galaxy sur- 
veys allow one to trace the evolution of stellar systems 
from high redshift to their present-day counterparts 



(e.g.lWolf et al.ll2004tlBorch et all2006tlFontana et al 



20MHcoville et al. 2006a). The combination of deep 



imaging surveys with the Hubble Space Telescope 
(HST) and ground based spectroscopic surveys play 
a special role, as they provide both structural infor- 
mation and stellar content for galaxies to early epochs 
dBrinchmann et alJll998UScoville et al]|2006bl) . 

Large area shallow surveys (e.g. SDSS) have 
provided a detailed census of galaxies in the local 
volu me by luminosity, stella r mass, and stellar con 



tent (Kauffmannet al. 2003; Nakamura e t all 12003 
Tremonti et alJl2004lMandelbaum et al.ll2006l) . There 



is a remarkably strong correlation between the struc- 
tural properties of galaxies, as quantified by the Hub- 
ble type for example, stellar content as revealed by 
integrated spectra, and stellar mass. Galaxies with 
stellar masses > 3 x 10 10 M Q have colors and spectra 
indicative of primarily passive evolution, while lower 



t Present address: Dept. of Astronomy & Astrophysics, University of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637 
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mass galaxies have blue colors and sp ectral features in- 
dicati ve of on-going star formation (Kauffman n et al 



2003). The red sequence galaxies are spheroid dom- 



inated, the blue sequence galaxies are primarily disk 
and irregular systems. Understanding the origin of 
the tight connection between spectral and morpholog- 
ical classes and their mass dependence is critical to a 
complete picture of galaxy formation and evolution. 

Intermediate depth surveys with HST combined 
with ground-base spectroscopy have shown that the 
Hubble sequence is largely in place at z ~ 1 (e.g. 



Lillv et al. 1995; Abraham et al. 1996; Brinchmann et al. 



1998UConselice et alJ2005l:lKajisawa & Yamadal2005l) . 
At higher redshifts, deep surveys (e.g. GOODS) reveal 
galaxies with structures th at are not easily classified 
into disks and spheriods dRavindranath et aL 2006; 
LotzetalJl2006l IPapovich et alj|2005b . This trend is 



even stronger in small area ultra-deep surveys e.g.the 



Hubb l e Deep and Ultra-Deep Fields) ( Abraham et al.. 
19961: lYan & Windhorst! 12004 ICoe et alJ l2006h . al- 



though some large disks were present a t z ~ 2.5 
dLabbe et al.l 12003: Sto ckton et ail |2004 . Red and 



Near-IR surveys have shown that the red sequence 
of massive galaxies was in place at z ~ 1, although 
with number dens ities lower than today by a fac 



tor of ~ 1 . 5-3 dBell et alj|2003t IChen et al.l 12002; 



Faber et al.l 120051) . Thus it appears that the galax- 



ies acquired their present day morphologies in the 
1 < z < 3 epoch. 

The onset of the strong correlation between mor- 
phology and stellar content and its dependance on 
mass is not as well constrained. One of the pri- 
mary goals of the G emini Deep Deep Survey (GDDS; 



Abraham et al.l 2004) was to use stellar-mass-selected 
samples to probe galaxy evolution in the critical 1 < 
z < 2 range. The GD DS and other samples (e.g. 



K20; iFontana et al.l 12004 have shown that the total 
stellar mass density evolves mildly for z < 1.5 and 
that the high-mass end in particular is slowly evolv 



ing dBrinchmann & Ellisll2000l; iDickinson et al.ll2003 . 
Bell et al.ll2004l:lGlazebrook et al.ll2004 . Many of the 



most massive galaxies at z < 2 have red optical- to 



near-I R colors dGlazebrook et al . 2004; Pap ovich et al 



2006 ) and spectra dominated by old stellar populations 
dCimatti et ap2004l: [McCarthy et aDl2004 . In this pa- 
per we examine the evolution of the stellar mass den- 
sity as a function of morphological and spectral types 
using the GDDS spectra and deep HST/ ACS imaging. 
Our analysis is aided by the use of quantitative mor- 
phological measurements obtained using a new 'quasi- 



Petrosian' threshold technique. This allows us to de- 
rive robust morphologies over a range of redshifts and 
apparent magnitudes. We show that the correlation be- 
tween morphology and spectral class remains strong 
for passively evolving objects at z ~ 1.5 and that 
the mass densities in spheroids (defined by morphol- 
ogy) and passive systems (defined by stellar content) 
evolve steeply in the 1 < z < 1.5 range. These re- 
sults strongly suggest that the core formation epoch for 
massive spheroids is drawing to a close in this interval. 
Signatures of major mergers in the passive population 
appear to have faded by z ~ 1.5 or higher, and subse- 
quent mass evolution at z < 1 apparently does not sig- 
nificantly disturb the spectral-morphology correlation. 
In a subsequent paper we will extend these results to 
z ~ 2 by using NICMOS to explore the rest-frame vis- 
ible morphologies of the most distant passively evolv- 
ing galaxies. 

We describe our GDDS+HST/ACS sample in Sec- 
tion 2, our approach to quantitative morphological 
classification in Section 3. Our results are presented in 
Sections 4 and 5, and discussed in Section 6. Our con- 
clusions are summarized in Section 7. An Appendix at 
the end of the paper presents a detailed analysis of the 
robustness of our morphological classification tech- 
nique. Throughout this paper we adopt a cosmology 
with #o=70 km/s/Mpc, tt M = 0.3, and il A = 0.7. 
The Vega magnitude system is used throughout. 

2. SAMPLE 

The GDDS is a spectroscopic survey of an optical 
and near-IR-selected sample targeting massive galax- 
ies at 0.8 < z < 2.0. The motivation for the survey, 
along with selection functions, sampling weights, de- 
tails of t he observations and c atalogs, are presented in 
Paper I dAbraham et al. I l2004 . Only a brief overview 
is given here. 

The GDDS sample is drawn from four fields in 
the o ne square-degree Las Campanas IR imag ing sur- 



vey (Mc Carthy et al.l 120011: IChen et all 120021) . span 



ning a total area of 121 square arcmin. Broad-band 
colors were used to pre-select galaxies likely to be in 
the 0.8 < z < 1.7 redshift range, and the spectro- 
scopic component of the survey emphasized obtaining 
high quality spectra for red galaxies at these redshifts. 
Very l ong exposures, using the 'Nod & Shuffl e' tech- 
nique dGlazebrook & Bland-Hawthornl 1200 lb to sig- 
nificantly improve subtraction of the night sky emis- 
sion lines, yielded high quality spectra from which 
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redshifts could be derived for 308 galaxies to a limit 
of Ivcga = 24.5 mag. The spectra define a one-in- 
two sparse sample of the reddest and most luminous 
galaxies near the I — K vs. I color-magnitude track 
mapped out by passively evolving galaxies in the red- 
shift interval 0.8 < z < 1.7. This sample is aug- 
mented by a one-in-seven sparse sample of the re- 
maining high-redshift galaxy population. The spec- 
tra go deep enough to allow redshifts to be obtained 
for L+ galaxies irrespective of star-formation history 
at z ~ 1.5, and the survey thus targets a mass-limited 
sample out to this redshift. Because of the importance 
of cosmic variance, the GDDS fields were carefully 
chosen to lie in regions of the sky where the number 
of red galaxies is near the global average (in order to 
avoid obvious over-densities and obvious voids). 

Forty five orbits of imaging observations with the 
Hubble Space Telescope's Advanced Camera for Sur- 
veys were obtained from August 2003 to June 2004. 
These orbits were distributed over 7 pointings with in- 
tegration times ranging from 11.7ks to 16.4ks (per 
pointing). In areas of overlap between pointings 
(which occur in our SA02 and SA22 fields), the ef- 
fective integration time can be much larger than this 
(32.6ks in SA02 and 36. Iks in SA22), albeit over a 
small area. Our rationale for obtaining relatively long 
integration times r elative to those obtaine d by surveys 
such as GOODS JGiavalisco etalJl2004 and GEMS 
( Rix et al. 20041) will be given in the next section. A 
tabular summary of our observations is presented in 
Table 1 . ACS images of starforming and post-starburst 
objects from this sample have al ready been presented 
in earlier papers in th is series (ISavaglio et al.l 2005; 
Le Borgne et al ] l2006b . while the sample as a whole is 
examined in the present paper. 

Our imaging data cover 67 square arcmin, or 55% 
of the total area of the GDDS, but by carefully choos- 
ing the locations of the pointings within the GDDS 
fields, and by varying the number of pointings per 
field, we were able to obtain ACS images for 63% 
of the galaxies in our spectroscopic sample. Fully re- 
duced ACS F814W images for each of the four fields 
are shown in Figure Q] Our field orientations were de- 
fined to produce additional depth for high priority tar- 
gets within the allocated orbit total. This resulted in an 
unusual and non-uniform coverage within each field. 
The objects within our GDDS+ACS sample are repre- 
sentative of the parent distribution in their local den- 
sities and distribution on the sky. Our sub-sample of 
GDDS galaxies imaged with ACS are neither drawn 



from unusually dense regions nor from low density 
voids. 

Using the methodology described in the next sec- 
tion, reliable morphologies can be obtained for galax- 
ies down to If814W = 24.0 mag. Our images con- 
tain 3596 galaxies in the LCIRS parent sample, 1519 
of which are brighter than our morphological analy- 
sis limit. Of these, 144 galaxies are in the GDDS. A 
detailed breakdown of the numbers of objects in indi- 
vidual fields is given in Table 2. 

3. METHODOLOGY 

The approach to quantitative morphological clas- 
sification adopted in the present paper is based on 
an updated version of the now fairly well-established 
technique of subdividing the galaxy population into 
classes on the basis of position on an asymmetry 
vs. co ncentration diag r am. This approach was first 
used in [Abraham et al. d 19961) to analyze the Hubble 
Deep Fiel d, and variations o f it hav e now been widely 
adopted ferinchmann et al~ 19981: Conselicd 2003; 



l Lotz etal.1 12004 IScarlata et al.l 120061 iLiskeretal 
120061) . While details of implementation differ amongst 
authors, the general approach has been shown to work 
well for subdivision into broad classifications (e.g. 
early-types vs. spirals vs. peculiars/mergers). Un- 
fortunately the more subtle distinctions between the 
morphological classes (e.g. elliptical galaxies vs. SO 
galaxies) are not captured by the system. 

Our analysis in this paper is based on a new 
publicly-available code, MORPHEUS, which produces 
a large number of measurements which can be used 
to constrain the morphological properties of galax- 
ies in addition to concentration and asymmetrjo The 
optimal set of parameters for characterizing galaxy 
morphology will be explored in Nair et al. 2007 (in 
preparation), but as will be shown below, for the sim- 
ple purpose of isolating early-type galaxies from all 
other galaxy types, a straight-forward diagram based 
on only two parameters, namely the Gini coefficient 
and asymmetry, works rather well. The Gini coef- 



1 M0RPHEUS can be downloaded from the following web site: 
http://odysseus.astro.utoronto.ca/~abraham/Morpheus Please note 
that the code base is evolving and the version used in the present 
paper is the July 2006 release. In addition to incorporating some rel- 
atively new parameters (e.g. the Gini coefficient from Abraham et 
al. 2003 and M20 from Lotz et al. 2004) MORPHEUS also incorpo- 
rates improvements suggested by others for ways to better measure 
well-established parameters such as asymmetry (e.g. Conselice et 
al. 2003). 
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Fig. 1. — Hubble Space Telescope Advanced Camera for Surveys images showing our SA02 field (top left), SA12 
field (top right), SA15 field (bottom left) and SA22 field (bottom right). As described in the text, the tiling of the field 
pointings were chosen to maximize the integration time on galaxies with known high redshifts. Details are given in 
Table 1. 
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Table 1 
Log of observations 



Pointing RA Dec # Visits # Orbits Integration time 

(J2000) (J2000) (s) 



SA02 (2 pointings, 14 orbits, 32676s total integration time, 19.9 arcmin 2 area) 

1 02:09:43.48 -04:36:42.0 3 7 16338 

2 02:09:37.43 -04:38:34.6 3 7 16338 



SA12 (1 pointing, 6 orbits, 14640s total integration time, 11.5 arcmin 2 area) 
1 12:05:26.71 -07:23:34.2 2 6 14640 



SA15 (1 pointing, 5 orbits, 12200s total integration time, 11.4 arcmin 2 area) 
1 15:23:50.38 -00:04:34.3 2 5 12200 



SA22 (3 pointings, 5 orbits, 36070s total integration time, 24.0 arcmin 2 area) 

1 22:17:38.15 +00:16:35.1 2 5 11670 

2 22:17:37.71 +00:14:21.6 2 5 12200 

3 22:17:45.09 +00:15:33.2 2 5 12200 



Table 2 

Summary of the ACS imaging sample 



Sample 


SA02 


SA12 


SA15 


SA22 


All 


Galaxies in the LCIRS parent sample 


998 


839 


748 


1011 


3596 


LCIRS galaxies with 18 < I < 24.0 mag 


375 


287 


315 


542 


1519 


Galaxies in the GDDS sample 


48 


29 


26 


90 


193 


GDDS galaxies with 18 < / < 24.0 mag 


39 


24 


11 


70 


144 
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fic ient is fairly new to a stronomy (being i ntroduced 
in Abraham et al. (2003) and improved in lLotz et al, 



(2004)), but it has been used in econometrics for nearly 
a century to quantify the ineq uality of w ealth distribu- 
tions in human populations ( Gini 1912b . When ap- 
plied to galaxy images, the Gini coefficient provides 
a quantitative measure of the inequality with which 
a galaxy's light is distributed amongst its constituent 
pixels, and it can be used as a sort of generalized con- 
centration index that does not depend on galaxy sym- 
metry. Furthermore, as shown the Appendix A, the 
Gini coefficient remains a surprisingly robust statis- 
tic even in the face of morphological A'-corrections. 
We therefore consider classification based on position 
in the Asymmetry-Gini diagram (referred to as the 
A-G diagram, or A-G plane) to be an updated and im- 
proved version of classification on t he basis of the 



Asym metry-Concentration diagram (lAbraham et al 
1996). The A-G plane will be the central diagnostic 
diagram in the present paper for purposes of morpho- 
logical classification. 



As noted earlier. lLotz et al.l (120041) presented a valu- 
able refinement of the original definition of the Gini 



coeff icient given in lAbraham et al.l (120031) . lLotz et al 



(2004) computed the Gini coefficient using circular 
apertures scaled to multiples of the Petrosian radius. In 
Section l3~2l below we will show how this idea can be 
extended using so-called 'quasi-Petrosian' isophotes, 
which retain all improvement s to the Gini coefficient 
presented by lLotz et al.l (120041) . and improve the statis- 
tic further so it works better for galaxies of arbitrary 
shape. But before describing our improved statistic, 
we will describe the importance of having data deep 
enough to be worth analyzing with these improved 
tools. Further discussion can be found in Appendix A. 

3.1. Required depth 

The long integration times obtained in our ACS 
fields were motivated by our need for reliable mor- 
phological classification of early-type galaxies out to 
z ~ 2. Morphological A'— corrections impact early- 
type galaxies in a qualitatively differe nt manner from 



how th ey impa ct late-type gala x ies (se e lAbraham et al 



19961) and Brinchm ann et al.1 dl998h for a more de 



tailed discussion). Late-type galaxies are dominated 
by irregular knots of star-formation which become 
more prominent as the observed wavelength probes 
further to the blue, typically leading to increased asym- 
metry. On the other hand, early-type galaxies are 
(by definition) smooth and rather symmetric, and the 



changes that occur are more uniform — the central 
parts of the galaxy remain brightest at all redshifts, and 
the asymmetry remains low. As is shown in greater 
detail in Appendix A, morphological classification of 
early-types from a single band of observation is pos- 
sible out to high-redshifts, provided the data go deep 
enough to probe the galaxies out to similar physical 
radii over the redshift range of interest. Obtaining data 
of sufficient depth is the crucial ingredient — prob- 
ably the greatest source of confusion in morphologi- 
cally classifying early-type galaxies at high redshifts is 
systematic misclassification of early-intermediate spi- 
rals as E/S0 systems, because at high-redshifts (1 + z) 4 
flux diminution lowers the surface brightness of disks 
below the threshold of visibility, thus increasing the 
prominence of the central bulge. Our ACS integration 
times were chosen to be deep enough to avoid this ef- 
fect. 

Figure |2] shows the signal-to-noise/pixel as a func- 
tion of age and redshift for HST ACS observations 
of a projected stellar mass density corresponding to 
that of the Milky Way at the solar radius. The panel 
on the left corresponds to an exposure time of 3ks, 
and the panel on the right corresponds to a 20ks in- 
tegration. Different lines correspond to observations 
of galaxies at z = 0.7, 1.2, and 1.7 (the range of red- 
shifts explored in the present paper). The logarithm of 
signal-to-noise/pixel is shown as a function of the age 
of the stellar population, assuming the fairly typical 
stellar population model described in the caption. To 
account (rather crudely) for the light smoothing and 
connected-component optimizations done by the pre- 
processing imaging segmentation software used (dis- 
cussed further below), we assume an effective pixel 
scale of 0. 1 arcsec/pixel, and that a SNR/pixel greater 
than 0.5 is detectabk^. On this basis, Figure [2] shows 
that at z ~ 0.7 one is able to probe out to large radii 
in the rest-frame of galaxies even with relatively shal- 
low (e.g. 1 orbit, or 3ks) optical observations obtained 
with HST. At this redshift, projected mass densities 
corresponding to the intermediate regions of large disk 
galaxies are observable for around 60% of the maxi- 
mum possible age for the galaxy (defined as the age of 
the Universe at the epoch of observation). However, 



2 This number is somewhat arbitrary, because in reality the limiting 
isophotal signal-to-noise is dependant on the smoothing parameter 
and connected-component algorithm (i.e. 4-point vs. 8-point con- 
nectivity) adopted by the segmentation software used to separate out 
a galaxy's pixels from the underlying sky. A value around 0.5 is rea- 
sonable assuming fairly heavy smoothing and 4-point connectivity. 
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3ksF814W 20ksF814W 




0123456 0123456 
Age (Gyr) Age (Gyr) 



Fig. 2. — Predicted F814W-band log 10 (signal-to-noise/pixel) as a function of age for a reference stellar population 
seen with the HST ACS camera at z = 0.7 (solid lines), z = 1.2 (dashed lines), and z = 1.7 (dotted lines). The 
panel on the left corresponds to an exposure time of 3ks, and the panel on the right corresponds to a 20ks integration. 
Calculations are based on a stellar population with a surface stellar mass density of 10 8 M Q /kpc 2 , which is comparable 
to that of the solar neighborhood. When applying K-corrections we assumed an exponential star-formation history 
with an e-folding timescale of 1 Gyr, a Miller-Scalo IMF, and 30% solar metallicity. To approximately account for 
the light smoothing and connected-component optimizations done by SExtractor when segmenting galaxy pixels from 
the background sky, we assumed an effective pixel scale of 0. 1 arcsec/pixel, and that a SNR/pixel greater than 0.5 
is detectable. (Note that log 10 (0.5) = —0.301). With these assumptions, the grey regions shown in the plot are 
unobservable due to low signal-to-noise. Another factor impacting the observability of the stellar populations is the 
age of the Universe, so the curves shown are limited to ages less than that of the Universe at each epoch of observation. 
See text for details. 
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beyond this redshift biases quickly become rather se- 
vere — structural information at optical wavelengths 
comes only from the densest regions in galaxies (such 
as bulges), or from the youngest parts of galaxies (such 
as bright star-formation complexes with ages of less 1 
Gyr). By z = 1.7 it is simply no longer possible to use 
shallow HST optical data to compare fairly fundamen- 
tal characteristics, such as the size and morphology of 
galaxies, with those of local objects. With 20ks inte- 
gration times, much fairer comparisons can be made 
— surface mass densities near that of the solar neigh- 
borhood with maximal ages can be probed in an un- 
biased way out to about z = 1.2, and at z = 1.7 
one is able to probe out to radii comparable to those 
probed by shallow HST imagery at z — 0.7. If, as 
suggested by Bouwens et al. 2005 (see also lye et al. 
2007), significant star-formation activity only begins 
around z ~ 6 (when the Universe is already about 1 
Gyr old, so that the maximal age of a stellar population 
at z — 1.7 is 2.8 Gyr instead of the 3.8 Gyr age of the 
Universe at that redshift), then Figure [2] suggests that 
one would be able to probe maximally old stellar pop- 
ulations at solar-neighborhood mass densities out the 
limits of our survey with 20ks integrations. Our goal 
in defining the rather complicated-looking ACS imag- 
ing field geometry shown in Figure Q] was to cover as 
many GDDS galaxies with known redshifts as possi- 
ble, subject to the constraint of having at least 12ks of 
integration everywhere, and up to 33ks of integration 
in areas where fields overlap. These areas of overlap 
were chosen to correspond to regions with many high- 
redshift galaxies. 

Having obtained data of the required depth, the next 
ingredient needed for reliable morphological classifi- 
cation of early-type galaxies is construction of a mor- 
phological catalog whose measurements have been 
synchronized (as closely as possible) to a common 
physical radius in the rest-frame of each galaxy. Stan- 
dard tools can take us part of the way to this goal, and 
the next stage in our analysis was to create a photomet- 
ric c atalog using SExtractor (IBertin & Arnouts 
Il996l) . which constructs a photometric catalog by iso- 
lating (or 'segmenting') galaxies from the background 
sky on the basis of fairly complicated criteria, based 
on differential thresholding followed by stages of im- 
age splitting based on connected-component analysis. 
Synchronizing measurements to a common physical 
rest-frame is done by using the SExtractor cata- 
log and segmentation maps as the input for our own 
software (MORPHEUS), which measures galaxy struc- 



tural parameters contained within a 'quasi-Petrosian' 
isophote, as defined in the next section. 

3.2. Definition of the Quasi-Petrosian Isophote 

The classic prescription for dealing with limiting 
isophote mismatches introduced by (1 + z) 4 surface 
brightness diminution of galaxies is to measure galaxy 
properties within a circular aperture who se size is 
a mu ltiple of a galaxy's Petrosian radius ( Petrosianl 
Il976l) . Our more modern formulation of this basic idea 
is valid for galaxies of arbitrary shape, and is thus more 
generally useful for analysis of the very diverse popu- 
lation of galaxies seen in deep HST images. 

The 'quasi-Petrosian' isophote is constructed us- 
ing an algorithm which works for galaxies of arbitrary 
shape. All pixels in the galaxy image (defined using 
the SExtractor segmentation map) are sorted in de- 
creasing order of flux to construct an array, fi, contain- 
ing the flux in the i th sorted pixel. This array is then 
summed over to construct a monotonically increasing 
curve of cumulative flux values: 



(i) 



In analogy with the definition of the Petrosian ra- 
dius for a circular aperture, we calculate the Petrosian 
isophote by determining the pixel index i which satis- 
fies the following equation: 



fi=V>< 



(2) 



Note that ^ is the cumulative mean surface bright- 
ness in the sorted array. The flux value of the pixel 
which solves the equation defines the quasi-Petrosian 
isophote. All pixels brighter than fi are retained, and 
all pixels fainter than fi are discarded. The parame- 
ter T] is fixed for all galaxies in the sample, and ex- 
perimentation has shown us that rj = 0.2 is a good 
choice, probing far enough into the outskirts of nearby 
galaxies to get outside of regions dominated by nuclear 
bulges, but not so far out that low and high-redshifts 
samples cannot be fairly compared. Because the data 
are discrete, Equation |2] will rarely be solved exactly, 
but a perfectly adequate approximation is obtained by 
simply noting the index of the first zero crossing of 
fi—rjx^^yA zero crossing is guaranteed to occur 
at the last pixel index, and when this occurs it indicates 
that no quasi-Petrosian isophote exists, and the algo- 
rithm converges on the standard isophote as the best 
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approximation to the quasi-Petrosian isophote. We 
consider this a rather graceful failure mode, and when 
it occurs, MORPHEUS simply flags this condition and 
proceeds with its analysis using the standard isophote. 

It is interesting to compare Equation[2]with the defi- 
nition of th e Petrosian radius a dopted by the SDSS col- 
laboration (iBlanton et al.l 120011) . In this formulation, 
one starts with Z(r), the azimuthally averaged surface 
brightness profile of the galaxy. For a given profile, the 



Petrosian radius, 



is the radius r at which which the 



following relation holds: 



/ o g 7 dr'2Trr'I(r')/{n(1.25 2 - 0.8 2 )r 2 
J r dr'27rr'/(r')/(7rr 2 ) 



0.2. 



(3) 

Thus the Petrosian radius is simply the radius at which 
the local surface brightness in a circular annulus is 
equal to 20% of the mean surface brightness within 
the annulus. Understanding this, we see that while 
Equation [2] looks very different from Equation [3] both 
operate in a similar way. The factor (50 is simply 
the mean surface brightness of a galaxy's brightest i 
pixels, so Equation|2]is essentially comparing the sur- 
face brightness at a point to a scaled surface bright- 
ness interior to an isophote — the key idea in the stan- 
dard definition of a Petrosian radius. If a galaxy's light 
distribution is circularly symmetric and monotonically 
decreasing from a central pixel, then the positions of 
the brightest i pixels will closely describe the filling-in 
of a bounding circle, so in this special case the anal- 
ogy between Equation [2] and the standard Petrosian 
formalism is nearly exact. However, Equation [2] re- 
mains valid for galaxies of arbitrary shape, and is both 
operationally easier to calculate and more robust than 
Equation[3] 

Figure[3]shows the the fraction of each galaxy's pix- 
els above the isophotal threshold defined by solving 
Equation|2]for r\ = 0.2 in our sample. As noted above, 
data that is too shallow for reliable classification and 
analysis is flagged by non-convergence of Equation|2] 
As shown in Figure [3] the fraction of each galaxy's 
pixels above the isophotal threshold is less than unity 
for all galaxies in our sample — therefore a quasi- 
Petrosian isophote exists for all galaxies in our sam- 
ple. This is not the case when our analysis includes 
galaxies fainter than Ifsi4W ~ 24 mag, and the ex- 
istence of a quasi-Petrosian isophote can be used as a 
sensible metric for determining how deep morpholog- 
ical measurements should be pushed to when analyz- 



ing a given datase{| It is interesting to consider how 
a sample defined by quasi-Petrosian area fraction re- 
sembles one based on a signal-to-noise thresholcQ, so 
the signal-to-noise vs. redshift distribution of the mas- 
sive galaxy sample considered in the present paper is 
shown in Figure [4] (See below for the precise defini- 
tion of this sample). As noted on the figure, our mas- 
sive galaxy sample has a mean signal-to-noise ratio of 
about 320. The subset of objects with z > 1.2 has a 
mean signal-to-noise ratio of about 220. 

Figure [3] also illustrates a drawback of the quasi- 
Petrosian approach. Because only a fraction of each 
galaxy's pixels are being used in our morphologi- 
cal computations, some information is being lost, al- 
beit the pixels with the lowest signal-to-noise in each 
galaxy. The median fraction of pixels lost is ^50% 
for galaxies at F814W < 22 mag, dropping to <20% 
for galaxies at F8 14W ~ 24 mag. Throwing data away 
is never a happy choice to have to make, but because 
of (1 + z) A cosmological surface-brightness dimming, 
it seems essential to do this (in a careful way) if fair 
comparisons of morphology over a range of redshifts 
are to be undertaken. Typically, only the high surface- 
brightness portions of nearer galaxies in a sample 
would remain visible at the highest redshifts probed by 
the same sample, so excising the low-surface bright- 
ness portions of bright objects is needed to harmonize 
regions of comparison over a broad range of redshifts. 

We conclude this section by noting, in passing, that 
there is a rather pleasing conceptual connection be- 
tween our definition of the quasi-Petrosian isophote 
given a bove, and the definitio n of the Gini coefficient 



given in lAbraham et alJ (120031) . Both statistics are ob- 



tained from simple computations on a ranked pixel list. 

4. THE MORPHOLOGICAL MIX OF GALAX- 
IES IN THE GDDS FIELDS 

Figure [5] shows small 'postage stamp' images of 
each GDDS galaxy, color coded by their spectral type, 
in the A-G plane. Galaxies with purely passive spec- 



3 In other words, if the quasi-Petrosian threshold is not being met for 
an object, then the data is so shallow that only the inner part of the 
galaxy is being probed. We thank Edward Taylor at Leiden Observa- 
tory for pointing out to us that meeting the quasi-Petrosian threshold 
condition is a convenient way to determine if data is deep enough 
for reliable morphological work. 

4 There are several different definitions for the signal-to-noise of an 
image. In the present paper a galaxy's signal-to-noise means the ra- 
tio of its background-subtracted flux to the product of the square root 
of its isophotal area times the standard deviation of the sky noise. 
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Fig. 3. — The quasi-Petrosian area fraction (the fraction of each galaxy's pixels above the isophotal threshold) plotted 
as a function /-band magnitude. 
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Fig. 4. — Signal-to-noise ratio as a function of redshift for galaxies with stellar masses log lo (M/M ) > 10.5. These 
objects constitute our massive galaxy subsample (see text for details). The mean signal-to-noise level of the galaxies 
shown, and of the subset of objects with z > 1.2, are noted near the top of the figure. 
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Fig. 5. — Asymmetry vs. Gini coefficient for galaxies in the Gemini Deep Deep Survey. The region of this diagram 
used to isolate early-type systems is marked at the bottom-right corner of the figure. Individual galaxies are shown as 
5 arcsec X 5 arcsec 'postage stamps'. Circled objects are included in the log 10 (M / M©) > 10.5 massive galaxy sub- 
sample defined in the text. The border of each postage stamp is is colored according to its spectral classification based 
on the system described in Paper I. Objects with spectral classifications corresponding to actively star-forming galaxies 
(classes 100, 1 10) are shown in blue; quiescent and nearly-quiescent systems (classes 001 and 011) are shown in red; 
intermediate activity systems (classes 010) are shown in green. Unclassifiable systems with known redshifts are shown 
with a gray border. (These unclassifiable objects had very limited wavelength coverage because they overlapped with 
nearby objects on our multi-object spectroscopy masks). Systems with no spectroscopic redshifts are shown without 
any border. See text for details. 
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tral types, identified with red borders, have primar- 
ily compact early-type morphologies and lie preferen- 
tially in the high Gini-low asymmetry region of Fig- 
ure [5] Galaxies with spectral types characteristic of 
active star formation have primarily disk and late-type 
morphologies and occupy a large range of the A-G 
plane, strongly favoring the high asymmetry and and 
low Gini regions. Galaxies with intermediate or com- 
posite spectral types occupy the region of the A-G 
plane between the passive and star forming galaxies. 
Figure [5] clearly shows that a strong correlation be- 
tween morphology, quantified by the A-G plane, and 
spectral class is in place at the redshifts spanned by the 
GDDS. For the early-type galaxies in particular, selec- 
tion by spectral class or morphology are nearly equiv- 
alent, and selection by morphology can be undertaken 
using simple linear cuts in the A-G plane. Through- 
out the remainder of this paper, we will define early- 
type galaxies to be those systems with G > 0.5 and 
A < 0.15. The choice of G > 0.5 is important, and its 
rationale is discussed in greater detail in Appendix A. 
The choice of A < 0.15 is based simply on visual in- 
spection of Figure[3] but is not particularly fundamen- 
tal for the limited purposes to which we are presently 
applying this diagram (namely, for determining the cu- 
mulative mass function of massive early-type galax- 
ies). This can be seen directly from Figure[5] where we 
have drawn circles around those objects which meet 
the mass completeness cut used when determining our 
mass functions (to be described in the next section). 
In fact, so few high mass systems lie above the asym- 
metry cut that our main conclusions do not change if 
no asymmetry cut is imposed at all. Therefore, our 
classification system for massive early-type galaxies is 
fundamentally based on the Gini coefficient. 

The concordance between morphology and spec- 
tral type is shown more clearly in Figure [6] where we 
present the individual images of galaxies classed as be- 
ing early-type using our formalism. Visual inspection 
of this figure is consistent with all galaxies being early- 
type systems. Approximately 15% (5/30) of the galax- 
ies do show weak disks, but without spiral structure, 
and these are probably SO or SO/a systems. Clearly the 
A-G diagram does not discriminate between E and SO 
galaxies. 



5. THE EVOLVING MASS-DENSITY FUNC- 
TION OF MASSIVE EARLY-TYPE GALAX- 
IES 

In this section we will attempt to synthesize the evo- 
lutionary history of massive early-type galaxies by ex- 
ploring the growth in their cumulative mass density in 
three broad redshift bins. The mass density functions 
were computed usin g the standard V m „, x form alism de- 
scribed in Paper III ( Glazebrook et al.1 2004), to whom 
the reader is referred for details. For present purposes 
it suffices to just highlight a few points that are impor- 
tant to bear in mind when interpreting the models. 

It is important to n ote th at our models assume a 
Baldry & Glazebrook! d2003l) Initial Mass Function 
(IMF), which has the same high-mass slope as a 
Salpeter IMF but which has a break at 1 solar mass 
(providing more realistic M/L values jfl Whenever 
we compare our results to those in the literature, we 
always convert to this standard IMF. Typical uncer- 
tainties in our model masses are around 0.2 dex in 
the main K < 20.6 sample considered in the present 
paper. 

Because the GDDS is a sparse-sampled survey, 
each galaxy in the survey acts as a proxy for a number 
of other objects with similar magnitudes and colors, 
so determining the stellar mass contributed by an in- 
dividual galaxy is only the first step in computing the 
mass function. The volume-weighted contribution of 
this galaxy to the integrated mass density in a redshift 
bin is computed using the sampling weights tabulated 
in Paper I. These weights account for the fact that the 
GDDS is not completely homogeneous, because when 
designing spectroscopic masks a higher priority was 
given to red galaxies than to blue galaxies. The sam- 
pling weights quantify the selection probability as a 
function of / — K and K, with reference to the full 
wide-area LCIRS tiles which formed the basis for the 
GDDS sample (and so also account approximately for 
the effect of cosmic variance). 

For maximum robustness in the present paper we 
focus our analysis on cumulative mass functions for 
systems above a mass threshold. We adopt a threshold 
of log 10 (M/M Q ) > 10.5. As shown in Figure 1 of 
Paper III, this is the limit to which the GDDS is mass- 
complete for z < 1.7. In other words, red galaxies 



5 The following convenient relationship (accurate to within a few per- 
cent independently of the SFH) can be used to convert between 
mass-to-light ratios measured using our chosen IMF and those de- 
termined with a Salpeter IMF: M/L K {BG) = 0.55M/ L K {SP). 
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Fig. 6. — Postage-stamp images showing the morphologies of the 30 objects classified as early-type galaxies on 
the basis of the system described in the text. Galaxies are shown sorted in order of decreasing redshift. Each im- 
age is 5 arcsec by 5 arcsec in size, and labeled with the galaxy's ID number, spectroscopic classification, redshift 
confidence class, rest-frame (U-B) color, redshift, and stellar mass inferred from our best-fit model. Objects without 
high-confidence spectroscopic redshifts have their redshifts labeled in parentheses. The borders of each image are 
colored according to spectroscopic classification, as described in the caption to Figure [5] Gray regions surrounding 
groups of postage stamps indicate which of three broad redshift bins the objects fall within. These bins are used to 
calculate the cumulative stellar mass function described in $5] 



SA1 5-7241 ClaSS=001 
(U- B) - 0.26 eonf-4 



SA1 2 6456 
(U-B)= 0.32 



class=00l 
conf=4 



SA22-0B96 
(U-B)= -0.01 



class=000 
confsO 







mass=JQ.700 z= 0.749 



mas5= 10.720 z=0.S12 



mass=-1.606 (7=0-530) 



14 



with stellar masses of \og 10 (M / M Q ) < 10.5 would 
begin to drop out of our survey at the highest redshifts. 
Another issue that arises when computing the cumula- 
tive mass function is the treatment of objects without 
redshift information. Because the GDDS has aredshift 
completeness of over 85 percent in the redshift range 
explored in this paper, the issue is not of overwhelm- 
ing importance, and none of the conclusions in this 
paper change if galaxies without redshifts are simply 
omitted. However, for the sake of consistency we will 
treat these objects in the same way as they were treated 
in Paper III. We therefore ass ign photometric r edshifts 
(computed for the LCIRS by lChen et all d2002ln to ob- 
jects without redshifts, and to objects whose redshift 
confidence class is < 2. (See Table 3 of Paper I for a 
detailed description of these confidence classes). 

The cumulative mass density function for systems 
more massive than log 10 (M/M Q ) = 10.5 is shown 
in Figure [7] The mass density function shown was 
computed in three redshift bins (0.8 < z < 1.1, 
1.1 < z < 1.4, and 1.4 < z < 1.7). The cor- 
responding error bars were estimated by a bootstrap 
analysis. Plot symbols in Figure [7] denote the follow- 
ing subsets of galaxies: (1) Spectroscopically quies- 
cent objects (spectral class 001) shown as diamonds; 
(2) morphologically selected (on the basis of position 
in the A-G plane) early-type galaxies shown as circles; 
and (3) all galaxies shown as squares. At the z = 
position we also plot the corresponding local stellar 
mass densities, bas ed on the analy sis of SDSS obser- 
vations given by iBell et al" ( 2003 ) and converted to 
our IMF, and incorporating our mass cut. (Note how- 
ever that the local early-type point is based on color 
selection, rather than on morphological classification; 
this point is discussed further below). T he cumulative 



mass f unct ions for all galaxy t ypes from lFontana et al 
(120041) and lBundv et al] ( 12006b are shown in gray. 



The most striking aspect of Figure |7]is that the cu- 
mulative mass locked up in morphologically-selected 
early-type galaxies appears to be increasing with cos- 
mic time substantially more quickly than the global 
rate. The cumulative stellar mass locked up in massive 
galaxies as a whole is increasing slowly by compari- 
son, changing by only about 0.3 dex over the redshift 
range probed, in agreement with the trends shown in 
Paper III. (This is reassuring, because our ACS imag- 
ing observations cover only 55% of the area in the 
complete GDDS). The cumulative stellar mass in spec- 
troscopically quiescent systems is consistent with be- 
ing fiat, although the error bars are large, and we can- 



not rule out that the spectroscopically quiescent sys- 
tems are evolving at a similar rate to that of the total 
galaxy population. 

In our highest redshift bin the fraction of the total 
stellar mass contained in massive early-type galaxies 
is similar to that contained by spectroscopically quies- 
cent systems (around 30% of the total). In our low- 
est redshift bin over twice as much mass is contained 
in early-type galaxies as is contained in quiescent sys- 
tems — in fact, early-type systems in our lowest red- 
shift bin contain about 80% of all the stellar mass 
locked up in galaxies with log 10 (M/Mq) > 10.5. 
This is consistent (within the considerable error bars) 
to the corresponding value in the local Universe (69%), 
which can be obtained from Figure [7] by simply di- 
viding the values shown for the two points at z = 0. 
Figure |7]therefore suggests that over the redshift range 
probed by the GDDS we are witnessing the formation 
of early-type galaxies as the dominant members of the 
high-mass end of the galaxy population, the same po- 
sition they maintain up to the present day. The result 
from the mass function is in entirely consistent with 
the visual impression one obtains by simply plotting 
the morphologies of massive galaxies in our sample 
as a function of redshift, as shown in Figure [8] The 
postage stamp images in this figure are shown in order 
of decreasing redshift, with early-type systems (based 
on our quantitative definition) circled. The steady in- 
crease with cosmic epoch in the fraction of early-type 
systems as a fraction of the massive galaxy population 
is rather striking. It seems fairly plausible to conclude 
that at at lowest redshifts probed by the GDDS, early- 
type galaxies comprise a similar fraction of the mas- 
sive galaxy population as is seen locally. 

6. DISCUSSION 

In the local Universe early-type galaxies are nearly 
always associated with quiescent stellar populations, 
so the steep evolution of early-types in the cumula- 
tive mass function shown in Figure |7]might be viewed 
as a surprise given the relatively shallow evolution 
in totally quiescent galaxies. A plausible explana- 
tion for this mismatch is offered by the possibility 
that morphological transformations and stellar mass 
assembly operate on different ti me scales. For ex- 
ample, it has long been argued (BauehetaD 1996: 



Kauffmann & CharloJl998l:lDe Lucia et al.ll2004 that 
the bulk of the stars in quiescent early-types galax- 
ies might have formed before these stars were orga- 
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Fig. 7. — Mass density functions for subsets of massive galaxies. As described in the text, our stellar mass cut is 
set at log 10 (Af/A7 Q ) > 10.5. Plot symbols denote the following subsets: quiescent (spectral class 001) galaxies 
[diamonds]; morphologically selected (on the basis of position in the A-G plane) early-type galaxies [circles]; and 
all galaxies [squares]. The c orresponding loc al stellar mass densities for massive galaxies, taken from the analysis of 
SDSS observations given by lBell et alJ (120031) and converted to our I MF and mass cut, are also shown. The two points 
shown in light gray are taken from DEEP2 survey work presented by lBundy et al. (2006), and from the K20 analysis 
published in lFontana et al.1 ( 120041) 
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Fig. 8. — Postage-stamp images showing the morphologies of the 54 galaxies in our sample with log 10 (M/Mq) > 
10.5, sorted in order of decreasing redshift. Early-type galaxies are circled. Each image is 5 arcsec X 5 arcsec in 
size, and labeled with the galaxy's ID number, spectroscopic classification, redshift confidence class, rest-frame (U- 
B) color, redshift, and stellar mass inferred from our best-fit model. Objects without high-confidence spectroscopic 
redshifts have their redshifts labeled in parentheses. The border of each galaxy image is colored according to the 
galaxy's spectroscopic classification, as described in the caption to Figure [5] Gray regions surrounding groups of 
postage stamps indicate which of three broad redshift bins the objects fall within. These bins are used to calculate the 
cumulative stellar mass function described in the text. 
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nized into spheroids, in which case the progenitors of 
such galaxies might be found in other morphological 
types. On the other hand, it might also be argued that 
the signatures of star-formation (e.g. blue colors and 
emission lines) remain in place for a period of time 
subsequent to the formation of an early-type galaxy, 
in which case the progenitors of quiescent early-type 
systems might still resemble early-types, but with blue 
colors and/or emission lines. We suspect that both sce- 
narios are at play in driving the steep evolution in the 
cumulative mass function of early-type galaxies shown 
in Figures [7] and [8] 

To further explore the link between morphology 
and star-formation history, Figure [9] shows the rest- 
frame (U — B) color of our sample, plotted as a func- 
tion of redshift. The rest-frame colors of our galaxies 
were computed using the best-fit spectral energy distri- 
bution templates used to construct our mass functions. 
The line shown near the top of the figure is the evo- 
lutionary track of an instantaneous starburst forming 
all of its stars at z = 3. Morphologically-classified 
early-type galaxies are circled. A number of inter- 
esting points emerge from this figure. As noted ear- 
lier, most morphologically-defined massive early-type 
galaxies are very red, but the figure shows that a num- 
ber of massive early-type galaxies with blue colors, 
probably from young stellar populations (as they do 
not exhibit AGN features in their spectra). The ex- 
istence of these objects should come as no surprise, 
since studies of internal color variations in early type 
galaxies have shown that a significant f raction of early- 
type galaxies at z ~ 1 have blue cores ( Abraham et alj 
1999; iMenanteau etaD 12004 IPapovich et al.1 120051) . 
Furthermore, such systems emerge from any scenario 
in which early-type galaxies are being built up over the 
redshift range of our survey with minor bursts of star 
formation occurring in time-scales longer than the stel- 
lar relaxation time. It is interesting to note that at high 
redshifts the red-end of the galaxy distribution lies 
consistently to the blue of the passive evolution model, 
and that at low redshifts most blue early-type systems 
are not very massive — presumably both these effects 



are m a nifestations of cosmic downsizing dCowie et al, 



19961: iKodama et all [20041 



Glazebrook 
Fontana et al.1 12004t iBauer et all 120051 



dCowie et al. 
et all 12004 
Juneau et al.1 



2005; iFaber et alJl2005UBundy et alJl2006h . We spec- 



ulate that most of the 'blue' early-type galaxies are ul- 
timately destined to turn into the quiescent early-type 
galaxies of intermediate mass. We also note that Fig- 
ures[5]and[8]both show that the abundance of this class 



of object would increase significantly if we were less 
stringent in requiring that early-type galaxies have low 
asymmetry, since inspection of these figures shows 
that a number of high Gini systems with high asym- 
metry have not been included in our sample of early- 
type galaxies. The asymmetric component of these 
massive objects is near their nuclei, and we suspect 
they are late-stage mergers. A number of such systems 
exhibit quiescent spectra, and these wou ld seem to be 
be go od candidates for "dry merging" ( Ivan Dokkum 
2005), though they may also be systems with high 
dust contamination. In any case, these objects are 
heavily outnumbered by star-forming systems. Fig- 
ure [9] also shows that a sizeable fraction of galaxies 
with very red rest-frame colors are not early-types. 
These objects comprise only 15% of the red galaxies 
at 0.7 < z < 1.0, but contamination rises quickly with 
redshift, and by z ~ 1.5 about half the objects with red 
rest-frame colors are not early-types. While the exis- 
tence of these objects is not a surprise — a number of 
recent papers have shown that color selection used to 
define early- type galaxy samples results in a mixed bag 



of ga laxies dYan & Thompsonll2003l : iMoustakas et al. 
120041) — the rapid change in the contamination rate as 
a function redshift is rather striking. 

Since color-selected samples cannot be used to de- 
fine samples of early-type galaxies without significant 
contamination, it is interesting to consider whether 
better samples of early-types might be constructed us- 
ing spectroscopy. FigurefTolshows the F814W images 
of the 21 galaxies classified as purely quiescent (class 
'001' in the taxonomy of Paper I) in our ACS sam- 
ple. About 2/3 of the sample (13/21) are early-types 
using our automated classifier. About half the remain- 
der are objects that would be classified as early-types 
if our asymmetry cut were slightly less rigorous, and 
on the whole the sample appears fairly homogeneous. 
We conclude that a fairly good sample of early-type 
galaxies could be defined simply by selecting purely 
quiescent systems. However, Figure [TT] makes it clear 
that care would have to be taken to exclude objects 
showing even small amounts of star-formation. This 
figure shows a montage of the 12 galaxies in our imag- 
ing sample with spectral classifications corresponding 
to a dominant old stellar population contaminated by 
low-level star formation (class '011' in the taxonomy 
of Paper I). It is obvious that these objects have a wider 
range of morphologies than those with the '001' (qui- 
escent) spectral class and a higher fraction of disks, 
and 2/3 of these objects lie outside the region in the 
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Fig. 9. — Rest-frame (U — B) color of our sample, plotted as a function of redshift. Galaxies with \og 10 (M / Mq) > 
10.5 are shown with large symbols. Morphologically classified early-type galaxies are circled. A fraction of the other 
red galaxies are massive are highly-reddened star-forming galaxies (Noll et al. 2006). The line shown near the top of 
the figure is the evolutionary track of a massive instantaneous starburst forming all of its stars at z = 3. See text for 
details. 
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Fig. 10. — 5 arcsec by 5 arcsec images of the 21 galaxies in our ACS imaging sample with spectral classifications 
corresponding to quiescent stellar populations (class 001). Early-type galaxies are circled. Gray regions surrounding 
groups of postage stamps indicate which of three broad redshift bins the objects fall within. These bins are used to 
calculate a cumulative mass function, as described in $5] See text for details. 
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Fig. 11. — 5 arcsec by 5 arcsec images showing the morphologies of the 12 galaxies in our imaging sample with 
spectral classifications corresponding to a dominant old stellar population augmented by low-level star formation 
(class 011). Early-type galaxies are circled. Gray regions surrounding groups of postage stamps indicate which of 
three broad redshift bins the objects fall within. These bins are used to calculate a cumulative mass function, as 
described in S|5] 
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A-G plane used to isolate early-type systems. Approx- 
imately 50% of the old + weak/truncated star forma- 
tion objects are disks and many show prominent HII 
regions. The remaining objects are compact and have 
early-type morphologies. We can only speculate on the 
ultimate destiny of the large disk galaxies in this sam- 
ple. Several of these objects have prominent bulges, 
and fading of their disks might bring them rather close 
to meeting our criterial for early-type galaxies. It is 
tempting to identify these with fairly recent merger 
events that will relax to the passive and early-type sys- 
tems at z < 1. Indeed these objects are closely re- 
lat ed to the massiv e pos t-starburst systems described 



Le Borgne et al. (2006) (Paper VI). They appear 



to have undergone significant star formation episodes 
at z > 1.5 and, barring further star burst episodes, 
will appear as passive systems on the red sequence by 
z ~ 0.8 — 1. (We note that spectro s copy o f color- 



selected samples - Yan & Thompson! (2003) - have 



also revealed luminous disk galaxies with red colors 
and spectra dominated by old stars.) 

Our results paint a picture of the morphological as- 
sembly of many massive early-type galaxies over the 
redshift range 1 < z < 2, but it is important not 
to take this conclusion too far. Large area surveys 
(Brown et al. 2006) clearly indicate that the most mas- 
sive systems have formation redshifts, both stellar and 
dynamical, well above the redshift limit of our present 
sample. The age analysis carried out in Paper IV and 



similar analyses of other samples ("e.g. |Cimatti et al. 



2004; Heavens et al. 2004; Jimenez et al 



2006) 



sug- 
gest that these objects formed at z > 3 on average and 
in some cases at considerably higher redshifts. As we 
will show in a subsequent paper, the rest-frame R-band 
morphologies of the passive galaxies with z > 1.3 
are consistent with early dynamical formation as well. 
These results are another manifestation of the popular 
down-sizing paradigm in which massive galaxies form 
early and less massive systems are ass embled at later 
times dCowie et al.lll996l: ljuneauet aDl2005h . At the 
low-redshift end of our sample, the convergence in the 
fractional mass density in massive early-type galaxies 
(relative to the total mass density at the same redshift) 
to something similar to the value in the SDSS indicates 
that the epoch of massive early-type formation may be 
drawing to a close by z ~ 1 (a conclusion consistent 
with the high post-starburst fraction described in Paper 
VI). 

The offset in the mass density contained in the low- 
est redshift GDDS point in Figure|7]and the local data 



point is rather striking, and this may indicate that what- 
ever process is driving the process of massive galaxy 
assembly at z < 1 is operating in a way which con- 
serves the overall fraction of early-types even as it 
builds the total stellar mass in the galaxy population 
from z = 0.7 to z = 0. However, we remain a little 
skeptical about this conclusion. The most robust state- 
ments that can be made from Figure|7]are based on rel- 
ative measurements internal to each dataset. Compar- 
ing the high-redshift GDDS points to their low-redshift 
counterparts in Figure [7] is dangerous, even though 
we chosen to p lot them together in the same diagram. 
Bell et al. (2003) note that there is an at least - 30% 



systematic uncertainty in the local estimate, and the 
uncertainty in the high-redshift points is likely to be at 
least this big. A nother source of potential bias is that 
Bell et al.l (12003b define their early-type sample on the 
basis of color, rather than on the basis of morphology. 
Much of the offset between the cumulative stellar mass 
locked in early type at z — and z > 0.7 may be due 
to systematic differences in the way samples of high- 
redshift and low-redshift early-type galaxies are con- 
structed — once again, we emphasize that the trends 
seen internally to either data set are more reliable. 

7. CONCLUSION 

We have described an improved methodology for 
morphologically classifying galaxies, based on 'quasi- 
Petrosian' image segmentation. This results in more 
robust measurements of galaxy properties at high red- 
shifts. This methodology has been applied to define a 
sample of early-type galaxies in HST ACS images of 
the Gemini Deep Deep Survey fields. Using this sam- 
ple, we computed the cumulative stellar mass func- 
tions of morphologically-segregated subsets of galax- 
ies with log 10 (M /Mq) > 10.5. Although our error 
bars are large, we find striking evidence for evolu- 
tion in the fraction of stellar mass locked up in mas- 
sive early-type galaxies over the redshift range 0.7 < 
z < 1.7. This redshift range corresponds to that over 
which massive early-type galaxies morphologically as- 
semble, and over which the strong color-morphology 
correlations seen in the local Universe begin to fall into 
place. 
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Appendices 



A. ROBUSTNESS OF THE QUASI-PETROSIAN GINI COEFFICIENT AS A GALAXY CLASSIFIER 

In this Appendix we will investigate the robustness of our galaxy classification methodology. As described in 
Section 4, our classifications are based on a galaxy's position in the Asymmetry vs. Gini diagram. For purposes of 
computing the mass function of early-type galaxies, we have shown in the text that the Gini coefficient plays a much 
more significant role than does asymmetry. In fact, the main conclusions of this paper would remain unchanged if 
asymmetry were neglected completely. Therefore our aim in this appendix is to demonstrate that the Gini coefficient 
measured in a single red band (F814W) can be used as a simple and very robust classifier of early-type galaxy mor- 
phology out to z ~ 2. We emphasize at the outset that this claim is true only when data is deep enough to allow 
quasi-Petrosian isophotes to be used to calculate the Gini coefficient, as described in the text. 

0.25 i ■ — ■ — ■ — ■ , — ■ — ■ — ■ — ■ — i — ■ — ■ — ■ ■ — , — ■ — ^ 
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Fig. 12. — Asymmetry vs. Gini coefficient for 800 nearby galaxies observed in g'-band as part of the Sloan Digital 
Sky Survey. Quasi-Petrosian thresholds were used when measuring both quantities, as described in the text. Galaxies 
classified visually as being early-type systems are shown in red. The vertical line denotes the G = 0.5 cut used to 
distinguish early-type galaxies from all other systems. 

Our analysis of the robustness of the quasi-Petrosian Gini coefficient is based on measurements of this quantity 
obtained through u\ g', r' and i' filters for 800 nearby galaxies in the Sloan Digital Sky Survey (SDSS). This is 
a subset of the morphological sample analyzed by Nair et al. (2007). We refer readers to that paper for additional 
details. Nair et al. (2007) will also present detailed comparisons of classifications made using quasi-Petrosian Gini 
coefficients to those made using other methods. Our ACS imaging sample has a signal-to-noise ratio floor of 100 (only 
a single galaxy in our ACS imaging sample has has a signal-to-noise ratio below this) so our analysis will be restricted 
to galaxies with a u'-band signal-to-noise ratio > 100. 
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Figure [T21 shows the Asymmetry vs. Gini coefficient diagram for g'-band imaging of our SDSS sample. Galaxies 
classified visually by one of us (Preethi Nair) as being early-type systems are shown in red. The vertical line denotes 
the G — 0.5 cut used to distinguish early-type galaxies from all other systems. It is seen that this simple bifurcation of 
the A-G plane into two parts does a rather good job of distinguishing early-type galaxies from all other systems. As 
expected, there is some leakage of visually-classified early-type galaxies to regions of the diagram with G < 0.5, and 
some leakage in the other direction (non-early-types into the G > 0.5). These systems are almost all early-type spirals 
whose visual classifications are known (based on comparisons of visual classifications between observers, e.g. Nairn 
et al. 1995) to be ambiguous (e.g. Sa vs. SO/a vs. SO systems). 

The G — 0.5 cut evidently does a good job of distinguishing high signal-to-noise g'-band images of SDSS early- 
type galaxies from the rest of the galaxy population, but how robust is this cut to changes in rest wavelength of 
observation, and to decreasing signal-to-noise? Figure [13] is an attempt to address this question. The top row of 
Figure [T3l shows the Gini coefficients of our SDSS reference sample measured at different wavelengths (using u\ g', 
and i' filters) plotted against the Gini coefficients in g' -band. The diagonal line shown in each panel delineates a perfect 
mapping between the two parameters, i.e. no change in Gini as a function of wavelength. The quasi-Petrosian Gini 
coefficient is seen to be remarkably robust to changes in rest wavelength when measured from data with reasonably 
high signal-to-noise. The Gini coefficient measured in g'-band barely changes when measured in 7-' -band and z'-band, 
and for early-type galaxies in particular the changes are hardly bigger than the random measurement errors (Nair et al. 
2007). Measurements of the u'-band Gini coefficient probe wavelengths blueward of the 4000A break, and our naive 
expectation was that Gini coefficients should be significantly different in u'-band when compared with measurements 
made at wavelengths redward of the break. At first inspection Figure Qj] does seem to show this, in the form of 
systematic offsets between the g'-band and it' -band images. However, a more careful inspection shows that the trends 
seen are almost certainly due to the the low signal-to-noise ratio of the u'-band data, instead of being due to systematic 
differences in the intrinsic g'-band and u'-band images of galaxies. 

Each panel in Figure [l3]records the mean signal-to-noise level of the sample being plotted. It is seen that the SDSS 
w'-band data has a mean signal-to-noise ratio of ~ 130, which is a factor two to three lower than the mean signal-to- 
noise ratios of our ACS images, and about a factor of five lower than the signal-to-noise ratios of the corresponding 
SDSS g', r', and V images. It is clear that these differing signal-to-noise ratios need to be accounted for before 
comparing the w'-band Gini coefficients of the data obtained at different wavelengths. This is a crucial point, because 
the ACS F8 14W filter begins to probe blueward of the 4000A break at z > 1, where the bulk of our data lies. Since the 
SDSS u'-band filter data would seem to provide the best match to the rest-frame wavelengths being probed by much 
of our ACS data, understanding the reason why the Gini coefficients of w'-band data shown in Figure [T3l appear be 
systematically lower than those measured at longer wavelengths is of special significance. The columnar base of the 
'T' in FigureQ~3]is an attempt to understand the w'-band data by showing how the Gini coefficient of g'-band sample 
changes when noise is added to the images to lower the mean signal-to-noise ratios of the samples plotted. The mean 
signal-to-noise ratio of the noise-degraded g'-band data in middle panel is similar to the mean signal-to-noise ratio 
for galaxies at z > 1.2 in our sample. The mean signal-to-noise ratio of the lowest panel is similar to that of the data 
in the SDSS u'-band sample. The strong resemblance between this panel and the u'-band panel at the top left of the 
Figure suggests that the apparent offset in the Gini coefficients of the high signal-to-noise g'-band data and the low 
signal-to-noise it' -band data is mostly due to poor signal in the latter, and not due to morphological K-corrections. 

A more explicit comparison of it' -band and noise-degraded g'-band Gini coefficients of our local calibration sample 
is shown in Figure [14] The mean signal-to-noise ratio of the g'-band sample has been harmonized to be nearly 
identical that that of the u'-band sample. Within the scatter introduced by random measurement error, Gini coefficient 
measurements at g'-band and u'-band are identical, even though these bands straddle the 4000A break (at which point 
visual morphologies of galaxies can appear to make large changes). This highlights the strongest benefit of using an 
abstract quantity like the quasi-Petrosian Gini coefficient as measure of morphology. This is also the essence of our 
contention, made repeatedly in the main text of this paper, that if quasi-Petrosian Gini coefficients are used to quantify 
morphology, then a basic 'cure' for the deleterious effects of morphological i^T-corrections is simply to obtain deeper 
data. 

We conclude this appendix with Figure Q3] which shows an object-by-object postage-stamp image montage of 
galaxies near the G = 0.5 border of the Gini coefficient vs. redshift diagram. This diagram can be used to inspect 
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Fig. 13. — Analysis of the robustness of the quasi-Petrosian Gini coefficient, based on measurements of this quantity 
obtained through u', g\ r' and V filters for 800 nearby galaxies in the Sloan Digital Sky Survey (SDSS). The x-axis 
of each panel in this figure is keyed to the g'-band Gini coefficient measured for each galaxy in this sample. The top 
row shows the Gini coefficients obtained at different wavelengths (through u , g\ and V filters) plotted against the 
Gini coefficients in g'-band. The diagonal line shown in each panel has a slope of unity and intersects the origin, and 
therefore corresponds to no change in Gini as a function of wavelength. Vertical and horizontal lines shown in each 
panel at Gini coefficients of 0.5 correspond to our proposed cutoff for early and late-type galaxies. Changes in the 
Gini coefficient are remarkably small over the entire wavelength range spanning the g' through i filter set. While 
Gini seems to change significantly when comparing u'-band to z'-band in the SDSS, this is easily shown to an artifact 
introduced by the low signal-to-noise in the SDSS u'-band data (which is typically much shallower than that obtained 
in other bands). To explore the sensitivity of the Gini coefficient to noise, the base of the "F in the figure shows how 
the Gini coefficient of g' -band sample changes when noise is added to the images. The mean signal-to-noise ratio of 
the noise-degraded 5' -band data in middle panel is similar to the mean signal-to-noise ratio for galaxies at z > 1.2 
in our sample. The mean signal-to-noise of the bottom panel is similar to that of galaxies with the poorest data in 
our high-redshift sample, and also to the mean signal-to-noise of the data in the SDSS it' -band sample. The Gini 
coefficients of noise-degraded g'-images and w'-band images are compared explicitly in the next figure. 
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Fig. 14. — Gini coefficients of u'-band galaxies in the SDSS compared to Gini coefficients of corresponding noise- 
degraded g' -band images. The mean signal-to-noise ratios of the samples have been harmonized to be nearly identical 
(~ 130). The scatter is consistent with the random measurement error at these low signal-to-noise levels (see Nair et al. 
2007 for details). As described in the text, when a quasi-Petrosian formalism is adopted for making the measurements, 
an individual galaxy's Gini coefficient is nearly independent of its wavelength of observation. 
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Fig. 15. — An object-by-object comparison of Gini coefficient vs. redshift for galaxies near the border of the G = 0.5 
cutoff threshold used to discriminate between early-type galaxies and all other systems. Objects in the gray region 
shown are worth looking at individually because they are at redshifts where F814W imaging is probing blueward of 
the 4000A break and also because they are near enough to the G = 0.5 cutoff that they could be brought below it at 
very low signal-to-noise. However, all these objects are seen to be high signal-to-noise systems, and only one object 
has a visual morphology even remotely consistent with that of an early-type galaxy. As in Fig. 5, the border of each 
postage stamp is is colored according to its spectral classification based on the system described in Paper I. See caption 
to Fig. 5 for details. 
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the morphologies of objects near the cutoff threshold used to discriminate between early-type galaxies and all other 
systems in the present paper. Objects in the gray region shown in Figure [15] shown are of particular interest, because 
these systems are at redshifts where F814W imaging is probing blueward of the 4000A break, and also near enough 
to the G = 0.5 cutoff that they could conceivably be brought below it at very low signal-to-noise. However, all these 
objects are seen to be high signal-to-noise systems, and only one object (with a high star-formation rate spectrum) has 
a visual morphology even remotely consistent with that of an early-type galaxy. 
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